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The author discusses particular solutions of a second order equation designated by source equation. 
This equation is special because the metric of the space where it is written is influenced by the 
solution, rendering the equation recursive. The recursion mechanism is established via a first order 
equation which bears some resemblance to Dirac equation. In this paper the author limits the 
discussion to solutions with constant norm but makes use of 4-dimensional hypercomplex numbers 
in matrix representation, a concept that is formally introduced in a section devoted to that aspect. 
The particular solutions that are found exhibit symmetries that can be assigned to spin, electric 
and color charges of elementary particles, leaving mass as a free parameter. Massless particles 
can also be assigned to special solutions of the source equation, with the cases of photons, gluons 
and gravitons clearly identified, together with another massless particle which does not seem to be 
related to anything detected experimentally. Another section deals with particle dynamics under 
fields, showing that both gravitational and electrodynamics can be modelled by geodesies of the 
spaces whose metric tensors result from the recursion mechanism. Finally the author suggests two 
lines of future work, one deriving fields from densities and currents of masses and charges and the 
other one aimed at determining particles' masses. 

PACS numbers: 02.10.Yn Matrix theory; 02. 20. -a Group theory; 12.10.-g Unified field theories and models. 



I. BACKGROUND 



One can conceive of a very large computer with the sole 
task of solving a small number of simultaneous differen- 
tial equations, with the only itch that the equations are 
evolutionary, meaning that in the course of time their 
form evolves as a result of previously found solutions. 
Assuming that the problem could be appropriately for- 
mulated, this would create an evolutionary system of 
equations whose complexity would surely outgrow the 
computing power of any machine. In a science-fiction 
pass one could also assume the computer to be itself a 
product of the equations, such that its computing power 
would grow in parallel to the equations' complexity. The 
author proposes that the Universe as a whole is such a 
system, whose evolution can be understood up to our 
ability to understand the basic equations. We will deal 
with evolutionary equations by writing non-linear equa- 
tions in 4-dimensional space; the solutions to these will 
be seen as evolutionary when 3-dimcnsional space is con- 
sidered. To this effect we will focus on one particular 
form of the basic equations as such equations are not 
known to us in general, although we can discuss some of 
their characteristics. This particular form is designated 
by source equation in Ref. [lj because it is the source of 
dynamics and dynamic space, an essential concept intro- 
duced in the reference cited earlier. We believe the source 
equation to be a condensation of the basic equations in a 
similar way to which Klein-Gordon and Hclmholtz equa- 
tions are condensations of Dirac and Maxwell equations 
0, ■ The reader may also find useful a look at Ref. Q 
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as this is a precursor of the present paper, establishing 
most of the assumptions used here, including setting the 
problem in a very specific 4-dimensional space, with the 
justification for such an approach. 

This discussion is about a mathematical problem, i. 
e. without any preexisting physical assumptions: All 
the variables in the equations are pure numbers, unless 
clearly stated otherwise, so the question of dimensional 
homogeneity is never an issue. This facet of the issues at 
hand is detailed in Ref. fj| where it was shown that di- 
mensions and physical interpretation can be introduced 
at a later stage, simultaneously with the definition of the 
physical constants c, G, h and e. 

The general form of the source equation is 

s^v^ - -i>, (i) 

where g^ v is a space metric, tjj is a function whose math- 
ematical characteristics will be discussed in the following 
section, and where V M designates the covariant deriva- 
tive with respect to the coordinate x* 1 . We use Einstein's 
summation convention, with Greek indices and super- 
scripts taking values between and 3 and Roman ones 
running from 1 to 3. It is presumed that the basic equa- 
tions of the dynamics will provide a means of determining 
g^ v from tp, in such a way that the successive equation- 
solution sets resulting from Eq. Q become coupled in 
a recursive manner, i. e. the solutions become self sus- 
tained. This process then links the space metric to the 
dynamics, defining what we have designated as dynamic 
space, a space which contains and is generated by the 
dynamics. The space metric will be considered separable 
in two components: an "inertia" reflecting the influence 
of the particular point on the dynamic space whose coor- 
dinates are under consideration, and a "field" reflecting 
the influence of all other points. 
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It is important to emphasize that dynamic space ex- 
cludes anything that does not participate in the dynamics 
under study. If we take the viewpoint of Physics for one 
moment, at the Universe's scale everything participates 
in the grand dynamics and is contained in the Universe's 
dynamic space. When we decide to study some particu- 
lar dynamics isolated from the rest of the Universe, this 
dynamics defines a space for itself, which excludes ev- 
erything else. Influences from other dynamic spaces can 
be accounted for via the field component of the metric. 
Consider the Earth's orbit around the Sun: We have the 
choice of using the dynamic space defined by the two bod- 
ies, or the dynamic space of the Earth alone, influenced 
then by a field component due to the Sun. These two 
approaches are equivalent if within the dynamic space 
concept we take universal space as a superposition of 
dynamic spaces, within which we can choose the com- 
ponents for the need of a particular study. 



II. 4-DIMENSIONAL COMPLEX NUMBERS 

A complex number can be represented in the exponen- 
tial form by 

V> = me ie (2) 

or equivalcntly in real-imaginary form as 

rj) = m(cos8 + isin#), (3) 

with m a positive number called the norm or modulus 
of the complex number. Either form is well suited to 
represent harmonic functions of one variable relative to 
another by taking just the real part. In the second form 
the variable ip as a harmonic function would be depen- 
dent on the variable 9 by the relation ip = mcos#, with 
the first form representing solutions 

V> = me ±ivx (4) 



In connection with dynamic space within the above 
view of superposed spaces, we define also an observer 
space, which is by definition obtained through a coor- 
dinate change that removes the inertia component from 
the metric of the dynamic space under consideration by 
the observer. Observer space can be then seen as the re- 
sult of a continuum of superposed dynamic spaces where 
each point can be mapped to a point of one particular 
dynamic space via inserting an inertia component in the 
dynamics equation of that space (as discussed in Ref. Q). 
In Physics, imagine one galaxy as dynamic space: When 
the galaxy's dynamics is studied in isolation its dynamic 
space extends to infinity. With the inertia component 
related to mass density, and with this density tending to 
zero with distance, the conversion to observer space con- 
fines the galaxy to a limited region of observer space. A 
similar link between dynamic and observer spaces can be 
applied to planetary systems, to the atom, and eventually 
to particles. 



This paper is the first in a series which will examine 
some mechanisms of dynamic space definition in a grad- 
ual process which has the potential to lead to universal 
space, starting here with "elementary" particles. Next 
section deals with the mathematical characteristics of the 
function tfj as the solution of Hclmholtz equation, also 
known as the "standing wave" equation, itself identified 
as contained in the source equation. We examine first 2- 
and 3- dimensional settings, with the ultimate intent to 
obtain a 4-dimcnsional formalism adequate for expressing 
the source equation solutions and metric. Once such so- 
lutions are found, their physical interpretation with and 
without fields follows to confirm their adequacy by lead- 
ing to key features of gravitation and electromagnctism. 
Particles involved with chromodynamics will only be in- 
troduced, without the necessary dynamic details. 



of Hclmholtz equation in one dimension 
d 2 

— $ = -v^ (5) 

with 9 = vx. 

A complex number is also appropriate for representing 
2-dimcnsional rotations since such rotations are harmonic 
relations between two orthogonal directions of space. In 
this interpretation a variation of coordinates in one of 
the directions is represented by the real part, while the 
coefficient of the imaginary part deals with the orthogo- 
nal direction. If the two coordinates are associated as a 
complex number 

x = x 1 + ix 2 7 (6) 

the rotated coordinates are given by 

x' = tpx, (7) 

which can be expressed in complex form as, 

x' = x 1 cos 9 - x 2 sin 9 + i (a; 1 sin 9 + x 2 cos 9) . (8) 

Complex numbers with norm unity belong to a group 
usually designated by £7(1), isomorphic to the group of 
2-dimcnsional rotations SO(2), whereby rotations in op- 
posite directions are represented by plus and minus signs 
added to the exponent of their exponential form as shown 
in Eq. J2J), and a 2n rotation in cither direction is equiv- 
alent to no rotation at all, i. e. the null rotation, with the 
exponential replaced by unity, and is commonly referred 
as a "full" rotation. 

In order to generalize complex numbers we designate 
by imaginary unit any traceless matrix u that verifies the 
condition 

u 2 = -I, (9) 
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then an imaginary matrix is defined as the product of a 
real number by an imaginary unit; the following relation 
is always verified 0: 



iu „tfu 



C = 



2 2 
I cos (9 + u sin 6*. 



(10) 



In a similar formal manner, we can define also traceless 
matrices h designated Hermitian units such that h 2 = I. 
For such matrices we have 



2 2 
= I cosh0 + h sinh0. 



(11) 



Rotations in 3 dimensions are associated with hy- 
percomplex numbers called quaternions (or 3D complex 
numbers) seen as an extension of complex numbers to 3- 
dimensional space [(|. We will define them through the 
Pauli matrices: 



1 

1 01 ' 



-i 

1 



, a 



1 

-1 



, (12) 



complemented with I as the identity matrix. They verify 
the well-known unitary and anti-commutator relations 



oV +oV = 2,1"' I. 



(13) 



with 5 13 being the Kronecker delta, making them Hermi- 
tian units as well as imaginary units when multiplied by 
i. 

When multiplied by i they do form a quaternion repre- 
sentational basis as defined by Hamilton 0] if the quater- 
nion product is simultaneously defined by: 



ia 3 <g> ia k 



a 3 a k , j f k, 
-I, j = k. 



(14) 



A unitary 2x2 matrix representing unitary 3D complex 
numbers has the form 



ip = c" 



(15) 



with the Oj 's real number components of a vector a with 
length 9 = (5 t3 aia^) 1 / 2 . From Eq. ((TUjl . Eq. lfT3|l is equiv- 
alent to 



(X ' (7^ 

xp = I cos + i-^- — sin 9, 



(16) 



which is a unitary quaternion, with the identity matrix I 
as the 0th order of its matrix representation. It is asso- 
ciated with a rotation of an angle 29 about the direction 
of unit vector a/0, as we shall see below. 

Defining the Hermitian coordinate matrix of a point 
(x 1 , x 2 , x 3 ) in 3D via the Pauli matricial basis 



i >J x 3 a 3 = I 



-x 2 + IX 1 



(17) 



we can see that the coordinates remain on a sphere of 
squared radius det 3 when transformed by the 

matricial similarity transformation 



(18) 



where ijr is the Hermitian conjugate of ip; notice that the 
second equality results from tp being a unitary matrix. 
The transformation above performs a generic rotation in 
3D containing the earlier 2D rotations with double their 
angle as we can see by taking a vector a along a 1 (only a\ 
is non-zero) and using Eqs. (|16[) and l|17|) in the similarity 
Eq. I|18(l . Coordinate x 1 is then unchanged, and the other 
coordinates transform as follows: 



x = x 2 cos(20) - x 3 sin(20), 
x' 3 = x 2 sin(20) + x 3 cos(20). 



(19) 



The same relation holds true when a circular permuta- 
tion of the indices is performed. This result means that a 
47r rotation angle is induced or two full rotations, for one 
cycle of the exponent 0, IE El ■ The group of 2 x 2 uni- 
tary matrices defined by Eq. (|15|) is designated by U(2), 
with SU (2) being the sub-group of those matrices whose 
determinant is unity 0, an d the group of real orthog- 
onal matrices of dimension 3, representing 3D rotations, 
is designated by 5*0(3). It is said that SU(2) is a double 
cover or two- fold covering of 5*0(3) in view of the angle 
doubling referred above jlfij . 

Generalization of complex numbers to 3 dimensions 
implies that we consider also non-unitary matrices; here 
the situation is more complicated than in 2 dimensions 
because the entities we are defining are 2x2 matrices 
and so must also be whatever generalizes the modulus of 
a standard complex. Eq. (|15(l can be generalized with 
the insertion of a traceless Hermitian matrix before the 
exponential; this factor is designated as modulus of the 
3D complex, which we will represent by \tp\: 

^ = \i'\c ia ^\ (20) 
The modulus relates to the norm H^H through 

H 2 = V^ = hA| t hA| = IIV'll 2 (2i) 



We will impose a further restriction on the modulus by 
taking the norm as a positive real number: 



(22) 



Now we are ready to verify that the 3D complexes de- 
fined through Eq. 1)2011 can be used as solutions of a 3- 
dimcnsional Hclmholtz equation. Consider the following 
equation with ip taken as a standard complex number 
function of coordinates x 3 and v a vector of components 



Vj and norm v: 



(23) 



usually the solution to this equation is written 

V = |V|e ±h ^\ (24) 
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where we have ignored an arbitrary phase factor and 



5 jk v-jV k . 



(25) 



This solution represents a plane standing wave with v 
seen as defining the orientation of the wave. We call 
standing wave to a standing harmonic 3-dimensional pat- 
tern. Imagine that ip represents light intensity: The 
space will then be seen as alternating between bright and 
dark along the wave direction; by defining planes normal 
to that direction spaced at half the wavelength we could 
split 3-dimensional space into alternating bright and dark 
zones. Using instead ifr as a 3D complex we can write dif- 
ferent solutions for the same equation 12311 : 



tp = |V|e ±h, ^< 



(26) 



We must verify that Eq. I|26|) together with Relations 
and (|2U is a solution of Eq. (|23J). The partial 
derivatives of ip are of the form 

djij) = ±|^f|iujtrj exp(±ii>fc<7 # ) = ±ivjajip (27) 

where we have lowered the index of o~j in order to avoid 
summation over j. The second derivative in x J of ip pro- 
duces a term <JjO~j, and considering Relations (|13f) . we 
get 



djjlp = - (Vj) xjj. 



(28) 



When the three second derivatives are added we obtain 
Condition JSHJ, and thus Eq. is verified. 

With either complex or 3D complex solutions for ip the 
general solution of Eq. l{2*3"|l is obtained through a super- 
position of positive and negative exponent solutions with 
different phase factors. The types of solutions, though, 
are different: Both equations represent standing plane 
waves but Eq. (|26|l is a associated to spinning of the 
points along the wave direction at twice the wave fre- 
quency; we will now call this a spin 1/2 wave. General 
solutions with 3D complex numbers are obtained by su- 
perposition of spin 1/2 waves; the norm of the 3D com- 
plex solution is interpreted as the wave amplitude, while 
the modulus contains information about its polarization. 
Spin 1/2 waves are especially indicated to describe axially 
guided waves, while no-spin waves, described by complex 
numbers, are better suited to describe unguided waves. 

The question remains about the possibility of defin- 
ing hypercomplex numbers representing solutions of 
Helmholtz equation in 4 dimensions. It is impossible to 
add a 4th imaginary unit on a par with the 3 quater- 
nion units we used and end up with a division algebra 
0; this does not mean, however, that we cannot define 
4-dimcnsional complex numbers. We will start with a 
discussion about 4-dimensional rotations, since in 2 and 
3 dimensions complex numbers have been associated to 
rotations and we are led to expect a similar relationship 
in 4D. 

Rotations in n-dimcnsional space are isomorphic to 
group SO(n) and determined by n(n — l)/2 parameters 



HUES, usuai ly chosen as rotation angles; this means 
that in 4 dimensions we have to select 6 different angles to 
characterize all possible rotations. There are 6 very spe- 
cial 4D rotations which preserve two of the coordinates 
while rotating the other two. For instance the orthogonal 
matrices inducing rotations that preserve x° , x 3 or x l , x 2 
are respectively 



(29) 




(30) 



When right multiplied by the column vector of the coor- 
dinates these special matrices rotate two coordinates by 
an angle while preserving the other two coordinates. 

Group SO {A) is doubly covered by group SU(2) x 
ST/ (2) in a similar way to the double covering of SO(3) by 
SU{2) 0; this suggests that we can replace Pauli a ma- 
trices by Dirac matrices built from them through matrix 
direct products. There are 15 of such matrices excluding 
the identity which can be generated by the formulas 



A J = I 2 0<7 J , 

p> = a j 0I 2 , 



(31) 



where I 2 is the 2x2 identity matrix and A B is the ma- 
trix direct product. Explicitly this set of Dirac matrices 
is given by 



A 1 = 



A 3 = 



P 2 = 



10 
10 
1 
10 



A = 



P = 



-i 

1 
-i 
i 

10 
1 
10 
10 







1 
i 



(1 
1 


Vo 







-1 

-1 



The complete set of 16 Dirac matrices is generated by 



(32) 



where X° = p° = I. 

The set of 16 Dirac matrices describes 15-dimensional 
space in the same way as the set of 4 Pauli matrices de- 
scribes 3-dimensional space. The source equation and the 
4D Helmholtz equation don't require this high dimension- 
ality because they can be associated to 4D rotations and 
hence to group SU(2) x SU (2), which can be represented 
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by Dirac matrices if appropriate relations are introduced 
between its elements, effectively reducing the dimension- 
ality from 15 to 4. The problem we will be facing is more 
complex, however because the source equation is obtain- 
able from a first order equation in much the same way 
as Klein-Gordon equation can be obtained from Dirac 
equation; the process is the same that was used in Rcf. 

From the cited reference we recover the first order 
equation, here in a slightly modified form 



(33) 

where are 4x4 Hermitian matrices defined by 

„o_ M (34) 



(35) 



where we define 4 matrices by 



a = E 10 , 



E 3j , 



which verify the relation 



25^. 



(36) 



Relation (|21() is used for the construction of s M . This is 
a more general procedure than was used in Ref. . 

Applying the operator is^V M to both members, and 
assuming the condition 



V^s" = 0, 



(37) 



the space metric of Eq. generated by the s M matrices 
is: 



-Ma" 



(38) 



The condition imposed by Eq. l|37f) automatically verifies 
Riemann's condition V ' a g^ v = although the converse is 
not necessarily true. 

There is an apparent inconsistency in Eq. (|38|l because 
s^s 1 ' 7^ s^s^ and thus the metric is defined as a non- 
symmetric tensor of 4 x 4 matrix elements. However, in 
many circumstances, as will be the case in this section, 
the metric will be anti-symmetric and cancellation be- 
tween corresponding terms will allow a diagonal metric 
to be used. 

Consider here expressions of the type g^V^, which 
imply an addition over all the \iv combinations; since the 
terms \w and v\x are symmetric and cancel each other, 
the expression is equivalent to what would be obtained 
with a diagonal metric. In index raising or lowering oper- 
ations it is not indifferent whether the metric is symmet- 
ric or anti-symmetric; we will always consider the metric 
to be defined by Eq. 138|) followed by elimination of anti- 
symmetric elements, making this elimination an integral 
part of dynamic space construction. 

In order to verify Eq. (|33|l and considering relations 
(|34|l we must try solutions of the type 



(39) 



where we have used ipo to denote the modulus and x° 
dependence of t\>. Solutions of the form given above can 
be inserted into Eq. I)33|l resulting in 



ilVlVoV* X>jV> 



(40) 



leaving only the x° dependence to be contended with. 

We will define 4D complex numbers to be 4 x 4 matri- 
ces, solutions of Eq. (|33|l . of the form 



ip = l^leP" 



(41) 



where up — la 3 and u° is some imaginary unit indepen- 
dent of the other 3. If we were only concerned with the 
simplified source equation above any Hermitian matrix 
with a trace of 4m could be used as modulus; this does 
not occur when we consider also Eq. J23J). For simplic- 
ity we start by examining the situations where po = m, 
Pj =0; equation l|33(l is simplified to 



is u V V' = ip- 



(42) 



Considering the definition of s° in Eq. (|34|) and noting 
that the covariant derivative of ip is equal to its partial 
derivative 



i| , 0|mu ?/' = m 2 ip, 



(43) 



where we have made | \ = ttl and multiplied both mem- 
bers by m 2 to remove the denominator. This equation is 
verified if 



\ip\ = imu°. 



(44) 



We are now in position to remove the simplifying con- 
dition inserted above and write the general solution 



ip = imu e 



(45) 



where s is a sign factor with the values ±1 and necessarily 
Condition l|52|l is verified. 

The extra imaginary unit u° needed to build a 4- 
dimensional complex number is built with recourse to 
matrices E 2 - 5 = p 2 X> . Any combination iqjp 2 \^ /n qi with 
qj equal to or 1 and (n q ) 2 the number of non-zero qj 
produces a suitable imaginary unit. Such combinations 
bear strong links with octonions [Tlj : Assuming justifi- 
able special rules on matrix multiplication that will en- 
sure non-associativity as required by the algebra, we can 
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define an octonion basis as follows 
ei = ip 2 X 1 , 

■ 2 \ 2 

e 3 = ip 2 A 3 , 

ip 2 (A 1 + A 2 ) 

V2 ' 
ip 2 (A 2 + A 3 ) 



e 4 



V2 ' 
ip 2 (A 1 + A 2 + A 3 ) 

ip 2 (A 1 + A 3 ) 
V2 ' 



(46) 



The multiplication table for octonions can be inferred 
from Fig. ^ where the arrows indicate directions of for- 
ward cycling, i.e. directions where sign is maintained; for 
instance e\e2 = e 4 or e 4 e6 = e^. In the reverse direction 
the sign of the result is negative. 



The justification for this rules is solely found on the fact 
that the elements of Eqs. (|46|l equipped with such prod- 
uct produce an octonion basis. As an example of use of 
the rules, we can check that e^e^ = —63: 

1. remove normalization factors, 

e 4 -» ip 2 (A x +A 2 ), 

e 6 - ip 2 (A 1 +A 2 + a 3 ), 

2. perform addition, 
ip 2 (2A 1 + 2A 2 + A 3 ), 

3. replace factors of 2 by zero, 
ip 2 A 3 -> e 3 , 

4. no rcnormalization needed, 

5. change sign for reverse direction — > — e^, 

6. not applicable. 

Under the assumptions above a 4D complex number 
has the form 




FIG. 1: Fano plane; the arrows indicate the directions of 
forward cycling; there are unseen arrows closing the cycles, 
for instance between and es or e$ and es. 

Using the multiplication table and the basis definition 
of Eq. H46f) we can establish the non-associative multipli- 
cation rule for octonions based on the a matrices: 

1. remove all normalization factors, 

2. perform the matrix addition of the two elements, 

3. replace any occurrence of 2A J by zero, 

4. renormalize to unit determinant according to the 
number of matrices in the sum, 

5. change sign if multiplying in the reverse direction 
as defined by the Fano plane, 

6. any clement squared follows matrix multiplication 
rules and equals —I. 



, PoQjP 2 ^ (\p q j p 2 X>x° k k 

ip = — cxp 1- rpfca x 

n„ V n„ 



(47) 



These numbers cannot be associated to transformations 
of 4-space because they involve all 15 dimensions. Pre- 
sumably the 15 coordinates are not independent and the 
actual space dimension is lower than 15 but this is a sub- 
ject which needs profound research which is beyond the 
scope of the present work. Some insight into the sort of 
transformations induced by the \a J imaginary units can 
be gained by considering the following matrix 



/ x° + ix 3 

? ■ 1 

x — ix 






-x 



x — ix 





ix 

3 



(48) 



ix 



Considering all the particular situation tp = exp(i0a 1 /2) 
and applying a similarity transformation coordinates a; , 
x 1 , a; 4 and x 5 are preserved while the other coordinates 



are transformed as 

a 

x - 

/3 

x = 

/6 

x - 
„/ 7 



x 2 cos ( 



-x 2 sin( 



x 3 sin ( 



I 3 cos ( 



x 6 cos — x 7 sin ( 



(49) 



x 6 sin ( 



x 7 cos ( 



If the other ia- 7 units were tested we would find the sets 
(x 1 , x 2 , x 3 ) and (a; 5 , x 6 , x 7 ) to undergo rotations of oppo- 
site sign. We are thus led to assume that coordinates x° 
to x 3 and x 4 to x 7 are interrelated, reducing the space 
dimension to 11; this is not an unexpected dimension and 
is coincident with what is used by superstring theory |l2j| . 
Transformations involving any of the u° units are more 
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complex, involving all the elements of the coordinate ma- 
trix. 

Considering Pj = in Eq. (|47|) , the unitary hypercom- 
plex numbers represented through the imaginary units of 
Eq. have the form 

where 9 = 2pox° is an arbitrary real number. If the 
upper left quarter of the coordinate matrix represents 4- 
space, as in l|48|) . these numbers represent 2 coordinate 
preserving transformations in the cases of e\, e2, e%, 1 
coordinate preserving transformations in the cases of e<4, 
es, e7 and no coordinate preserving transformations in 
the case of e§. These numbers can be classified according 
to the signs of p and qj , and according to the number of 
non-zero qj . For each sign of p there is one hypercomplex 
number for (n q ) 2 equal to 3 and 3 hypercomplex numbers 
for (n q ) 2 either 1 or 2. 

There is another unexplored possibility for the Oth 
imaginary unit which consist of making u° = ia°; this 
generates yet a different sort of transformation analo- 
gous to a global harmonic pulsation of the whole 4-space. 
Consider ip = exp(ia°d/2 and apply a similarity trans- 
formation to the coordinate matrix; all the elements in 
the upper left quarter, those representing 4-spacc, will 
appear multiplied by cos 9. 

These generalized complex numbers will be used in 
the present and following sections to write solutions of 
the first order and source equations. We may have used 
some discretion in choosing the mathematical entities we 
defined as 4D complex numbers; we were guided in this 
choice by need to make a connection with elementary 
particles later on. 

The source equation Q can become extremely intri- 
cate; this is to be expected, since it is applicable to all 
the problems of dynamics. The exposition of this sec- 
tion will be made clearer by a preliminary exploration 
of the 4D Hclmholtz equation's solutions using the 4D 
complexes introduced above. Consider Eq. I|23|) with j 
and k indices replaced by /x and v in order to make it 
4-dimensional 

S^d^ = -mfy (51) 

Just as in the 3D case, this equation accepts complex 
solutions of the type ?p = \ip\ exp(±ip A1 x AI ), with the con- 
dition 

S^PuPv = m 2 - (52) 

We are interested in exploring 4D complex solutions of 
the form given by Eq. I|41|) ; no matter which choice is 
made for u°, it will always be possible to write a solu- 
tion in the form of Eq. (|45() with Condition l|52|l verified. 
This solution represents a spin 1/2 wave in the direction 
of the unit 4- vector p whose components arc p^/m; m 
is also the wave frequency along the direction p; in 4D 



we will associate the designation Compton frequency to 
the wave frequency. The 4D wave is made up of two 
components: One 3-space component generated by the 3 
elements Pj, with frequency ^2{Pj) 2 , corresponding to the 
De Broglie frequency, and one component in the direction 
x° , with frequency po , corresponding to a transformation 
that does not preserve x°. When u° = ia° the x° wave 
induces a pulsation of 4-space, as referred. 

In this section we want to explore the simplest possible 
forms of the source equation, together with its solutions. 
The form of the equation is governed by the metric which 
includes an inertial component and a field component; we 
will assume here the latter to be the identity matrix, so 
we will limit ourselves to identifying the recursion mech- 
anism which determines the inertia component. 

If we consider 4D complexes defined by Eq. I|41|) . a 
subsequent application of Eqs. i|34|) and l|38|) generates 
an anti-symmetric metric given by the following relations 

oo, MM i 
g IMP IMP' 
m = M^ 
9 IMP' 

j0 a^M (53) 

9 IMP' 



y "IMP' 

9 jk = -9 kj =a j a k \W\ 2 

If | ip\ anti-commutes with a 3 this is a situation where 
the metric is anti-symmetric and can be replaced by a 
diagonal metric; making the substitution |M| ~~ ¥ m2 
one gets 

g» v = -^6' w , (54) 
m 

as the metric to be used in the source equation. 

Leaving aside the anti-symmetry conditions, having as- 
sumed the norm to be constant, and noting that a con- 
stant metric zeroes the Christoffcl symbols allowing the 
covariant derivatives to be replaced by partial ones, the 
source equation becomes the 4D Helmholtz equation 
(|5 1 p . Apart from a phase factor, this equation allows so- 
lutions of the type given by Eq. (|41J) with Condition l|52l) 
met and with diagonal modulus. If we want not only the 
source equation but also the first order equation (|33|l to 
be verified, then the modulus is given by Eq. (|44() . 

III. A PHYSICAL EVALUATION OF THE 
MASSIVE PARTICLE SOLUTIONS 

What is specially interesting about 4-dimensional 
standing waves in comparison to their 3D counterparts 
is that there are different sorts, according to the imag- 
inary unit u° that is chosen. Looking at the octonion 
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basis (|46fl . elements (ei,e2,e3) identify transformations 
where 2 of the 3-space coordinates are preserved, ele- 
ments (ci, es, ey) refer to transformations preserving only 
one 3-space coordinate and element ee 4D combined 
transformations involving all 3-space coordinates. For 
the 7 possible u° choices that form an octonion basis the 
concept of octonion units can be extended to the waves 
they generate, suggesting that octonion products per- 
formed among 4D waves may describe particles involved 
in chromodynamics. The different possible choices for u , 
including ia , point to the interpretation that Eq. I|45|) 
represents 8 particles and their respective anti-particles, 
with this solution describing the known elementary par- 
ticle dynamics. While the equations have so far been 
written as pure numbers, dimensions and physical in- 
terpretation are possible through the use of dimensional 
factors introduced in Ref. 0. 

Returning to Eq. 1)450. assume that pi = P2 = so 
that the particle moves along the x 3 direction. We have 
m 2 = (po) 2 + (753) 2 and the first thing we notice is that 
the frequency (p ) of the x° wave is reduced as (^3) in- 
creases, i. e. the transformation involving a; with 1, 2 or 
3 spatial coordinates, is transferred to rotation around 
x 3 . Secondly we notice that the Dc Broglie frequency is 
P3 with spin frequency of twice that value; we refer to 
this fact saying that the particle has spin 1/2; actually 
the spin can also be —1/2 when the sign factor s is nega- 
tive. Rewriting Eq. (|45(l for a stationary particle, defined 
by pj — 0, whose u° is chosen among the octonion basis 
elements 



, mqjp \ J ( simqjp \° 

tp = — exp - 

n q \ n q 



(55) 



we identify the particle's mass as m, connecting it with 
the norm of the particle's wave function; for a stationary 
particle the mass coincides with the frequency of the x° 
wave. 

The particle's electric charge can be defined as ^2qj/3 
and comes from the transformation involving the 0th co- 
ordinate through one octonion unit. For this reason, 
when (n q ) 2 equals 1 or 3, we will choose negative rather 
than positive qj, thereby allowing the electric charge to 
be —1 when n 2 — 3. Such a particle can be then associ- 
ated with an electron. 

The electric charge is +2/3 when n 2 = 2 and —1/3 
when n 2 = ljCorresponding to the up and down quarks, 
respectively @, El ■ Color charge is associated with the 
various qj 7^ possibilities and we can see that electric 
charge appears as a spherically symmetric color charge. 
Anti-particles are in turn represented by an equation 
identical to Eq. (|55|l where the sign of the qj's is reversed. 
The interpretation made above shows that Eq. (|55|l has 
all the necessary features to represent all leptons and 
quarks with the exception of neutrinos, with the parti- 
cle's mass being the only free parameter. It is of course 
desirable to remove all free parameters, deriving mass 
from some characteristics of the source equation's solu- 
tions, but we will not address the subject in this paper. 



We are going to associate neutrinos with the u° = ia° 
situation; as for the other particles, the neutrinos' masses 
appear as free parameters. 



IV. SOLUTIONS WITH ZERO NORM - THE 
QUESTION OF MASSLESS PARTICLES 

It is important to investigate solutions of Eq. (|ST|) when 
in = noting that the equation can be rearranged as 



Whenever 



<P*V,- fc ^ = -VooV'- 



VooV" = w i> 



(56) 



(57) 



and covariant derivatives can be replaced by partial ones, 
the equation above becomes the 3D Helmholtz equation 
(|23|l . for which there arc plane wave solutions of the type 
given by Eq. Q24[l. i. c. no-spin waves. Being in 4D, in- 
tegration of Eq. I|57(l has to be made with tp as a 4 x 4 
matrix; after integration we get 



^ = ^i,2,3e 



(58) 



where h is an Hcrmitian unit and "01.2,3 is a 4 x 4 matrix 
function of the 3 spatial coordinates. The general form 
of a plane wave solution for Eq. (|56|l can be written 



tp = Ce ±hulx0 e isp ' xJ : 



(59) 



with C an Hermitian matrix integration constant. 

Physical interpretation of Eq. I|59|) reveals a plane wave 
restricted to 3D space and evanescent in the 0th dimen- 
sion, if the positive sign is rejected. The phenomenon of 
evanescent waves is well known from total reflection and 
waveguides, where both exponent signs are also present 
but the positive sign is ignored for lack of physical sig- 
nificance. Here we interpret the solutions given by Eq. 
(|59|l as representing massless particles, of which we dis- 
tinguish three different types according to the Hermitian 
unit h that is present: An Hcrmitian unit obtained from 
the octonion basis element e$ by product with i carries 
electric charge in the evanescent wave and will be asso- 
ciated to photons; an Hermitian unit obtained from one 
of the other octonion basis elements by product with i 
carries color charge in the evanescent wave and will be 
associated to gluons; finally the Hermitian unit a carries 
no charge and will be discussed later in this section. 

Some words are needed here about the reason why we 
chose no-spin wave solutions for Eq. (|56[) while spin 1/2 
solutions exist. We could argue that spin 1/2 solutions 
cannot coexist in the ground state due to the Pauli exclu- 
sion principle, but we can do better than that: A no-spin 
wave defines a direction in 4-space but not a worldlinc; 
this is true both for dynamic or observer space. On the 
contrary, the observer space appearance of a spin 1/2 
wave is similar to a 4-dimcnsional waveguide and it de- 
fines a worldlinc coinciding with the waveguide axis; this 
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idea was first suggested in Ref . jj| . The number of cycles 
along any waveguide length is actually a measure of time 
interval in the waveguide's (or particle's) own clock. We 
can also assume all clocks in all waveguides to be syn- 
chronized from a common "zero time" in the origin of 
the Universe. When two waveguides meet in the same 
3D position and their time measurement coincides they 
have to interact; Ref. 01 includes the discussion of a 
collision problem in 4-dimensional optics, which is an il- 
lustration of this fact. The same does not happen with 
no-spin waves which never actually meet in one position 
because they are not associated to lines and this matches 
the behavior of massless particles. The interpretation of 
time as the count of cycles along the waveguide is en- 
tirely compatible with the use of dt as line element that 
we make in the next section. 

A superposition of an even number of ±1/2 spin waves 
has many characteristics of a no-spin wave; consider for 
instance the superposition of 2 equal modulus equal fre- 
quency +1/2 spin waves with different phase, both spin- 
ning around the x 1 direction: 



spins and multiplying the two wave functions we get 



i 1 ^ 1 _|_ ^ia. 1 (ujx 1 -j-p) 



(60) 



it is sufficient to make p = u for the equation above to de- 
generate into a no-spin wave given by ip = \ip\ cos(wx 1 )I; 
naturally the same could be done if the two waves were 
counter-spinning or with any even number of waves co- 
or counter-spinning. The spin of any wave will thus be 
expressed by the algebraic sum of the superposed waves 
spins. 

All solutions verifying Eq. I|59(l have complex depen- 
dence on the spatial coordinates, not quaternionic-like so- 
lutions. As we said before, complex numbers describe no- 
spin waves, whereby these solutions all exhibit no spin. 
Only combinations of an even number of particles with 
spin 1/2 can then lead to integral spin waves, behaving 
like no-spin ones. We shall examine this case next. 

We know from experiment that photons are associated 
to the annihilation of an electron and a positron; since 
we have established the wavcfunctions of the interven- 
ing particles it is elucidating to investigate how the an- 
nihilation process can be understood in terms of those 
wavefunctions. The wavefunctions of an electron and a 
positron are given by 



/ ^ — smeQX -\-'\r>~a? x 3 



smeQX +ipfca x 



(61) 



where ee is the octonion basis element corresponding to 
electron and positron and the exponent sign designates 
each particle's electric charge; the sign of s stands for the 
two possible spin orientations. By adding an imeex term 
within both exponents to represent the coupling between 



<g> V>+ => TO 2 e - S me 6 z (l-i) ^me^ (l+i) 



m 2 e 2s m ie 6 x- I 2 ^ e i Pj a- 



me 



'( 

ie6X ° (2 * v' p " ' j . 



2*e i; ' " '' 



(62) 

62T) 



In the first stage of the process (|62|) . the two particles 
are assumed to interact by means of the ia; terms in 
the exponents; in the second stage (|62l ) an intermediate 
particle is formed, which absorbs the result from the an- 
nihilation in the pj ; the superposition of the two waves is 
indicated by 2*. Finally the third stage 1)621' ) produces 
2 identical photons of spin 1 and frequency equal to the 
mass of one particle. The final result is a solution of 
Eq. (|56|l . representing a superposition of 3D spin waves, 
evanescent in the x° direction. As before only a negative 
evanescent exponent is acceptable and since s can have 
both signs the sense of the mutual orbiting must change 
according to the spin sign. 

The annihilation of two neutrinos, if at all possible, 
leads to a similar expression, where ie§ is replaced by a ; 
the resulting particles are 



ip = mc 



2*e'^ Q 



(63) 



This would be a massless particle not involved in any of 
the four known fundamental interactions and we do not 
dare proposing an interpretation for it, although it might 
be related to dark matter in the Universe. 

Returning to Eq. I|62|l . further annihilation is possible 
through the interaction of two i/'-V'-i- pairs with oppos- 
ing evanescent exponent signs s. We rejected the positive 
sign for its lack of physical significance, however two pairs 
with opposite signs produce pure 3-dimcnsional waves, 
not accompanied by any x° evanescence; these waves are 
the superposition of 4 wavefunctions and thus have spin 
2. Four identical particles with mass m can be annihi- 
lated producing 4 resulting particles represented by 



1> 



1 •: (•''" " ' 



(64) 



We believe these are good candidates for gravitons 
which may have been detected experimentally long ago 
by Allais' measurements of Foucault's pendulum anoma- 
lies and already discussed by the author 



SOLUTIONS WITH FIELDS 



In sections ITTT1 and HVI we discussed field free solutions 
of the source equation in situations which reduced that 
equation to the 4D Helmholtz equation; those solutions 
where found to have characteristics which suggested their 
physical interpretation as massive and massless parti- 
cles. However the real test of this interpretation relies 
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on the ability of the particles thus defined to react prop- 
erly to applied fields, accurately predicting the dynamics 
of those particles. Fields are assumed to alter the par- 
ticle's dynamic space through the s M matrices defined in 
Eq. (|34[) . The exact way in which fields are created is not 
yet clearly understood but we believe it is connected to 
the flux of massless particles; we think that the flux of 
masslcss particles determines the field gradient, the net 
value of the flux being linked to the norm of the field gra- 
dient and the flux direction to its direction. Things get 
more complicated when massless particles are originated 
in moving frames because an extra rotational factor in- 
tervenes. The whole process will eventually be governed 
by equations similar to Maxwell's equations albeit gener- 
alized to encompass the three field types: gravitational, 
electromagnetic and color. In the following exposition 
fields are presented axiomatically and the dynamics they 
induce are then seen to represent dynamics under the 
influence of gravity and electromagnetism; a loose asso- 
ciation is made with the massless particles identified in 
the previous section but no attempt is made to propose 
a formal link between particle flux and field gradient, 
which we intend to do in a future paper. 

We are going to assume that fields leave the particle 
wavefunction ip virtually unaltered, so that we can use 
the wavefunctions derived previously; this means that we 
can make = imu° and \ \ip\\ = to. It will frequently be 
convenient to write the modulus in the equivalent form 
IV^I = im cxp(±u°7r/2); as a first generalization step the 
s M matrices are modified from Eq. (|34[) by introducing 



G M , scalars, and A^, matrices: 



n iG M c^e ±u A W 2 
s = — - , 



G„ 



(65) 



Eq. I|33|) maintains its validity and by applying to both 
members the operator is M V M we get the source equation 
(Q) with g* v defined by Eq. J3B|) as before. 

In a general situation the metric defined by Eq. 138(1 
can be expanded to 



±£(^) 2 C ±U ° A ^ 



= 


1 




Tfl 


9 j0 = 


i 

9 




to 


9 ]k = 


1 

o 




n r 



(66) 



This metric will now have to be examined for various 
possible forms of G M and A M and the respective dynam- 
ics related to known interactions. Simultaneously a con- 
nection will be made to the massless particles identified 
in the previous section in order to associate them with 
mediators of the various interactions. 

The simplest case to examine is of course the one where 
all the are zero and the G M coefficients are all equal, 
G, t = G/4; inserting into Eq. I|66|) . the metric becomes 



anti-symmetric and we can reduce it to the diagonal met- 
ric 



5^. 



(67) 



This metric defines a space whose geodesies reproduce 
the dynamics of a particle with mass to in a gravitational 
field G, as long as the field sources are stationary. This 
problem was analyzed in Ref. and will be reproduced 
below for completeness. 

Inserting Eq. (|67[l into Eq. (Q, with if; given by Eq. 
(|41|l and assuming G to be a slowly varying function so 
that the Christoffcl symbols can be neglected 



G 



6^ Pfl p„ = 1: 



(68) 



then replacing^ = g^ v x u , with i M = drr^/di the- deriva- 
tive of the coordinate with respect to the geodesic line 
element, here designated by dt, 



/my 

Kg) 



1. 



(69) 



Here we remind the reader that although we have chosen 
the letter t for the line element, this is in no way associ- 
ated with time at this stage; the procedure whereby this 
association can be made later on was detailed in Ref. 0] . 

The equation above is called a dynamic space equation 
because the inertia element to is present in the metric; it 
can be converted into an observer space equation by the 



coordinate change X^ = mx^ . 

^ V X»X V 
G 2 



We can then write 



= 1. 



(70) 



The procedure to derive the 4 geodesic equations from 
Eq. (|7U|I can be found in many textbooks, for instance 
|17j . We define a constant Lagrangian with the value 1/2 
and insert into the second member of Eq. I|7U|I which be- 
comes 2L. Following that, we derive the Euler-Lagrange 
equations in a standard way. If G is independent from 
coordinate X° there is conservation of the corresponding 
conjugate momentum 



G 2 



1 

■ 

7 



(71) 



with 7 an arbitrary constant. For the other 3 coordinates 
we can write 




8 llv X> l X v djG 
G 3 

d G 

--^T, using J70J). 



(72) 



The reader is now referred to Ref. 0] in order to verify 
that in the particular case of a single body of mass M 
as origin of the field introduced in Eq. I|65|) . one must 
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make G = exp(— M/r), with r the distance between the 
two masses, thereby giving a clue that InG is the New- 
tonian gravitational field. In the cited reference it was 
also demonstrated that predictions of general relativity 
are verified to the first order approximation by the equa- 
tions above. 

The link of gravitational interaction introduced in Eq. 
(I65[l with the graviton equation (|64[1 can be made by 
noting that both G and graviton's modulus are products 
of the identity matrix by a scalar; it is thus arguable that 
the gravitational field can be mediated by such particles. 

Next we want to investigate the adequacy of function 
t/> to represent the electrodynamics of a charged particle. 
We start by taking G as unity and 



(73) 



with scalar; we will also assume that u° = 

±e§, effectively limiting the analysis to the elec- 
tron/positron family members. As a result the expo- 
nentials exp(±u°A M 7r/2) in Eq. (|66fl can be replaced by 
cxp(±A M 7r/2). It is also convenient to define 



(74) 



The metric elements become 



and V-x = x 1 qB 3 \ B 3 represents the magnetic field aligned 
along a; 3 and q = ±1. The Lagrangian is now defined by 



or o 

— = 8^x v - 2qx 1 B 3 x°x 2 + (qx^zx 2 ) . (77) 

777 V ' 

Some straightforward calculations lead to the conju- 
gate momenta 



Po 
m 2 
Pi 

9 

777 
P2 

m 2 

P3 

o 



x° - qx 1 B 3 x 2 



x\ 



x 3 . 



qx 1 B 3 x <) + (qx l B 3 ) 



2 .2 

X , 



(78) 
(79) 
(80) 
(81) 



Since the Lagrangian is independent from x 3 , Eq. I|81|) 

gives 



x = constant; 



(82) 



then, independence of the Lagrangian from x° and Eq. 
GHJ) give 



1 

7717 



qx 1 B 3 x 2 , 



(83) 



jo 



a j a k 



(75) 



This is no longer an anti-symmetric metric such as we 
have found in all previous cases because g 0j and g 3 do 
not cancel each other completely. After cancellation of 
the anti-symmetric parts the metric becomes 



5 kX V k V x Vi V 2 V 3 

Vi 1 ■ ■ 

v 2 • 1 • 

v 3 • • 1 



(76) 



The lower subscript metric g^ v is obtained, as usual, by 
calculating the inverse of g^ v 



(Vb) 2 



/ 


1 


-Vi 


-v 2 


-Vs 


\ 




-Vi 


(v ) 2 + (y!) 2 


V1V2 


V1V3 




-v 2 


V1V2 


(Vb) 2 + (V 2 ) 2 


V 2 V 3 




V 


-v 3 


V1V3 


V 2 V 3 


(V ) 2 + (V 3 ) 2 


1 



We don't develop here the details of electrodynamics gen- 
erated by the metric above; the procedure is similar to 
what was used for gravity and can be found in Ref. 0, 
where both electric and magnetic field dynamics were ex- 
amined. 

Consider now the situation where d fl V u =0 except for 
diV2 = qB 3 ; for simplicity we make Vq — 1, V\ = V3 = 



with 7 a constant greater than 1, unity applying to a 
stationary particle. 

Inserting Eq. l|g3*|) into Eq. i|SU|) and deriving the second 
member with respect to the line element At and equating 
to zero because the Lagrangian is independent from x 2 



qB 3 .! 

X . 

7777 



(84) 



Finally, from Eq. I|79|) , deriving the Lagrangian with re- 
spect to a; 1 and making the needed substitutions we get 



qB_ 

7777 



3 -2 
X . 



(85) 



Equations (|84() and l|85() represent the Lorentz force ex- 
erted by the magnetic field B 3 over a particle of electric 
charge q and moving with velocity x J . The frequency 
qB 3 /(m"f) is the angular frequency of the movement. 

In the limit 7 = 1 and x 3 = the particle approach 
breaks down because this approach is similar to the ge- 
ometric or ray approach of optics and cannot be used 
when any features of the worldline are of the same order 
of magnitude as the particle's Compton wavelength. If we 
take q — —1, for the electron case, the angular frequency 
becomes uj = —B 3 /(m'j) and one possible solution for 
Eqs. iJHUl and (JSSJl is 



x 1 = r cos (ut) . 
x 2 = r sin (u>t) . 



(86) 
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The use of t as an independent variable is here linked 
to a specific worldline in dynamic space; t can always 
be expressed in terms of the coordinates by inversion of 
the equations above t = &rct&ri(x 2 / x 1 ) / uj . Inserting Eqs. 
into Eqs. (|T5j) to l(ST|) we get the conjugate momenta 



m 

Pa = — , 

7 



Pi 



mB^r sin (u>t) 



1 



(87) 



P2 =0, 

P3 = arbitrary constant. 

Choosing zero for p 3 the particle's wavefunction can be 
written 

/ m n n mB 3 r sin (ait) . i , 

V> = exp — u u x u H — - u/x 1 

V 7 7 



exp 



u°a; + ia L m z ujr z sm z (wt) 



1„,2, 2 



7 



which is now an expression not involving x 3 and not 
breaking down when 7 — ► 1. In the limit there is spin 
precession with the cyclotron frequency o> = —B^/m; for 
a positron the frequency would show a positive sign. 

Similarly to the relation of gravitons with the gravi- 
tational field, photons can be seen as related to the EM 
field by their modulus when comparing Eq. I)73J1 with Eq. 
H62(l . In both cases there is an exponent which can in- 
teract with the exponent in the wavefunction modulus; 
photons can be considered as mediators of this field. Nat- 
urally relations like Eq. i|73[l can be established for the 
other octonion basis elements, originating the color fields 
and chromodynamics; this is a subject for a forthcoming 
paper. 

We cannot end this analysis without touching upon 
the situation where the G M differ from each other; this 
must surely be related to gravity, since we have seen how 
gravitational field is modelled by equal coefficients. We 
presume at this point that different G M 's must be associ- 
ated with mass currents in a similar way to what relates 
the Vfj, 's to electrical currents. This would mean that dif- 
ferent Gn coefficients model situations of non-stationary 
gravitational sources. As we said at the beginning of 
this section, it is believed that a set of coupled differ- 
ential equations links the components of each field and 
that these equations are a generalization of Maxwell's 
equations; we will investigate this possibility also in a 
forthcoming paper. 



in this work to the more traditional hyperbolic space- 
time. The character of their dynamical symmetry has 
been found to follow that of elementary particles, leav- 
ing mass as the only free parameter. Electric and color 
charge result directly from rotational symmetry proper- 
ties in 4D of standing waves, and massive particles were 
shown to verify the laws of dynamics following geodesies 
of a suitably defined 4-space. The physical origin of spin 
has been identified without appeal to Relativity princi- 
ples, and found to result merely from Helmholtz equa- 
tion requirements in 4D. The modifications to the source 
equation needed to express the complete dynamics were 
loosely tied to the character of massless particles, this 
being an area of further development as detailed below. 

Two complementary lines of future work are envisaged, 
following paths only suggested in the present work. The 
first line will try to complete the set of recursive equa- 
tions, deriving fields from particle densities and currents 
in dynamic and observer space, and in particular deriv- 
ing chromodynamics. As a mere working hypothesis we 
believe that a composite field can be defined as 



(summation is suppressed) (89) 



where G ;i represents the contribution from gravitational 
sources both stationary and non-stationary and A^ rep- 
resents both electromagnetic and color fields, as discussed 
in the previous section. The composite field will then be 
related to the flux of massless particles by a composite 
field tensor defined by 



(90) 



which is expected to verify relations similar to Maxwell's 
equations, namely 



8^F — T 



0. 



(91) 
(92) 



These are untested hypothesis and must be seen only as 
a working program. 

The second line of work will be a search for variable 
norm solutions of the source equation. Preliminary work 
has already shown that allowance for a variable norm 
and the consequent consideration of Christoffcl symbols 
can send the source equation into some modified forms of 
Emden-Fowler equation [T^ | which can determine mass, 
thus eliminating the only remaining free parameter. 



VI. CONCLUSION AND FURTHER WORK 
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